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Abstract 

We discuss systems containing a heavy quark and a heavy antiquark 
in the infinite mass hmit of QCD. Studying the limit of equal velocities 
for both heavy particles, we formulate an effective theory approach to 
heavy quarkonia-like systems. The method is well suited to processes 
in which the two heavy quarks annihilate, such as electromagnetic and 
strong decays of charmonium and bottomonium and weak decays of 
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1 Introduction 



In the past five years considerable progress has been made in heavy quark 
physics by studying systems with a single heavy quark in the infinite mass 
limit of QCD The mass m of the heavy quark sets a scale large compared 
to the intrinsic scale Aqcd of the light QCD-degrees of freedom, and the 
appearance of such disparate mass scales allows us to use an effective theory 
treatment of systems with a single heavy quark. This effective theory, the so- 
called heavy quark effective theory (HQET) has the interesting property of 
additional symmetries which are not present in full QCD. These symmetries, 
in addition to the usual machinery of effective theory, are a powerful tool 
in heavy hadron physics, which allows us to make QCD based and in many 
cases even model independent statements. The progress in this field is well 
documented in more or less extensive review articles 0. 

Almost all the applications considered so far deal with the one heavy 
particle sector of QCD. In HQET particle and antiparticle number are sep- 
arately conserved and all applications deal with either a single heavy quark 
or a single heavy antiquark. First attempts to consider states containing two 
heavy (anti) quarks or a quark and an antiquark revealed some difficulties 
when one calculates QCD radiative corrections. The anomalous dimensions 
of operators having matrix elements with states containing two heavy quarks 
turn out to be complex, at least if interpreted in the naive way. In addition, 
the imaginary parts behave as 1/v, where v is the relative three velocity of 
the heavy quarks. Subsequent investigations @] showed that the imaginary 
parts of the anomalous dimensions yield phase factors which have to be inter- 
preted as the non-abelian analogue of the Coulomb phase well known from 
electrodynamics. These phases are an infrared contribution, which has to 
be absorbed into the states. After redefinition of the states the anomalous 
dimensions are real, well behaved as the relative velocity becomes small, and 
hence are the true short distance contribution. 

Heavy quarkonia have to show up in the sector of HQET containing a 
heavy quark and a heavy antiquark. The velocities of the two heavy particles 
in such a state differ only by a small amount of the order Kqcd/^-, and 
hence one wants to describe heavy quarkonia in the limit, where the two 
heavy quarks move with the same velocity, which is then identified with 
the velocity of the heavy quarkonium. However, such a limit is ill defined 
for static quarks due to the divergent phases mentioned above. A proper 



1 



treatment of these phases hes at the heart of the formulation of a Heavy 
Quarkonium Effective Theory (HQQET), since these phases are related to 
the binding mechanism of the heavy quarkonium states. We shall see in what 
follows that HQQET may not simply be related to the two particle sector 
of HQET, because the divergent phases prevent us from taking the strict 
infinite mass limit. 

In this paper we shall clarify the underlying assumptions necessary to 
formulate HQQET based on the static limit. In the next section we shall give 
a theoretical description of the method, which has been applied to inclusive 
heavy quarkonia decays already in and compare to related ideas of Bodwin, 
Braaten and Lepage (BBL) 0]. In section 3 we discuss how to calculate QCD 
and higher order non-perturbative corrections in HQQET. 

2 Formulation of HQQET 

Our starting point is the Lagrangian of QCD, which we expand in inverse 
powers of 1/m 0, . The part of the Lagrangian involving the heavy quark 
is unique up to terms of order 1/m and is given by 

£ = U^\ivD)h[^^ + hi+\pP_t]phi+^ + 0(l/m2), (1) 

\2m/ 

while the corresponding expansion for the field Q of full QCD is given by 

Q = exp(-imvx) (l + — + 0(l/m^)] h'^+\x), (2) 
V 2m / 

where P± = (1 ± |^)/2. The superscript (+) indicates that the field 
describes a static quark moving with the velocity v; correspondingly we in- 
troduce also the field which describes an antiquark moving with the 
velocity w. The Lagrangian and the expansion of the full QCD field for the 
antiquark field is obtained from (|I]) and (0) by the replacement v —>■ —w. 

Let us first consider only the static, mass independent term of the ex- 
pansions (|1|) and (|^) and write the Lagrangian for a two-particle system 
consisting of a static quark and a static antiquark as 

£ = hi'-\tvD)hi^^ - hl;\iwD)hl;^ . (3) 
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Based on such a Lagrangian we may now consider the matrix elements in- 
volving two heavy particle states; as an example we shall study 



G= {A\Q{x)TQ{x)\0), (4) 

where A is a state containing a heavy quark and a heavy antiquark moving 
with velocities v and w respectively. In the static limit this matrix element 

becomes 

<^static = {A\hi"'\x)Thl~\x)\0)eMMv + w)x], (5) 

where the tilde denotes the static limit of the state. Logarithmic dependences 
on the heavy quark mass may be calculated using renormalization group 
improved perturbation theory in the framework of HQET. The one-loop QCD 
radiative corrections to a current of this type have been calculated in 0] and 
[^]. It turns out that in a naive calculation the anomalous dimension seem 
to acquire an imaginary part, which for v —>■ w develops a divergence of the 
general structure 

Im 7 = f{as) r — -= = (6) 
[vwY — 1 



where is known up to two loops 

/(«.) = (l + ^ 



31 10 ■ 



+ ■■■ • (7) 



The real part of the anomalous dimension vanishes as f — > w due to current 
conservation, and the solution of the renormalization group equation with a 
purely imaginary anomalous dimension yields for t; — > w a phase factor of 
the form 

1 r. /(«) ^^ ^ 



exp(i(/)(f«;)) = exp \ i f da-^ \ ; /?(«(//)) = 

(8) 

which is ill behaved in the limit v ^ w. 

This divergence of the imaginary part prevents us from taking the limit 
f — > w for two heavy static particles. On the other hand, this is exactly 
the limit in which we want to describe heavy quarkonia states. From this 
we conclude that the purely static Lagrangian (H) in not appropriate for the 
description of heavy quarkonia states. 
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In fact, the phases are the non-abehan counterpart of the Coulomb phases 
well known in QED. They are related to the long range part of the one gluon 
(one photon) exchange potential, which decreases too slowly and thus leads 
to infrared problems. Consequently, these phases are an infrared effect and 
have to be absorbed into the states [Q. 

In the channels, where this potential is attractive, bound states may 
occur, and these phases are directly related to the binding mechanism. This 
may be explicitly seen for the abelian case using eikonal methods, which 
correspond to the heavy mass limit [|l^ . The binding of a heavy quarkonium 



is clearly an infrared effect and has to be reproduced by the dynamics of a 
properly constructed effective theory for quarkonia. In such a state the two 
velocities differ only by a small amount of order 1/m which is a hint that we 
need to go beyond the static limit to describe quarkonia states. 

In order to see how higher-order terms of the Lagrangian cure the prob- 
lem, we make use of the fact that the full QCD Lagrangian is independent 



of the arbitrarily chosen velocity vectors v and w [^. The only kinematic 
quantity entering in full QCD are the true momenta of the particle p = mv+k 
and the antiparticle p' = mw — k', where k and k' are the residual momenta 
of the two heavy quarks. Thus we are led to define 

W iD 
V = f H and W = w -\ (9) 

m m 

corresponding to p/m and p'/m respectively. These combinations are invari- 
ant under infinitesimal reparametrizations of the velocities v v + 6v and 



w — > w+6w [^, since under such a reparametrization we have D — > D—m6v 
for the quark and D ^ D — m6w for the antiquark moving with the veloc- 
ity —w. In a full QCD calculation the singularity corresponding to the one 
appearing in (^ occurs in the imaginary part of the vertex function as 



Im T{p,p') = f(as)^====\n - , (10) 

where A now is an infrared regulator, revealing the infrared origin of the 
singularity. The vertex function of full QCD depends on the full quark mo- 
menta p and p', which are split into a large piece mv (mw) and a residual 
part k {k') as we switch to the effective theory. This is, however, an artificial 
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procedure, and the singularity may be reproduced either in the dependence 
on the velocity or through the residual momenta. 

In heavy quarkonia the velocities of its heavy constituents are almost 
equal, vw ~ 1, and one would rather go to the limit v = w and reproduce 
the divergence of the imaginary part as a singularity in the dependence on the 
residual momenta. This is achieved formally by reinserting the full momenta 
for the velocities in the divergent phase, i.e. by the replacement w — >■ V and 
w W, and we obtain from (Q) 

exp{i(t){vw))U+'>Th^_:2 — > /i(+)exp(#(V>V))r/iL~J . (11) 

If we now consider the limit v ^ w we have also V W, but now the phase 
depends on the residual momenta, which are represented by the covariant 
derivatives acting on the heavy quark fields 

(12) 

which means that we have rewritten the singular phases in such a way that 
they now depend on the difference of the residual momenta rather than on 
the difference of the velocities. 

However, the difference between v and V is a term of higher order in 
1/m, which was added in such a way that V is a reparametrization invariant 
quantity. In order to construct a leading order Lagrangian capable of repro- 
ducing the (infrared) phase factors and also generate binding of the two heavy 
quarks, we rewrite the static Lagrangian (^ in a reparametrization-invariant 
form; in this way we obtain as the leading order Lagrangian 

where we have replaced w —w in the last step. In the Lagrangian (|1^) 
we now may perform the limit v —>■ w without encountering a problem in 
the calculation of QCD radiative corrections; we now have already to leading 
order the scale of the heavy mass m, and ultraviolet contributions show 
up as ln(A/m) (A being the ultraviolet cut-off) and infrared contributions 
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exp(i(/)(VW)) exp(i(/)(V^)) = exp < 



da 



a., (m) 



(3a 



as ln(m/A) (A being the infrared cut-off). In this way a clean separation 
between (calculable) short distance effects and (non-perturbative) infrared 
contributions is achieved. The diverging phases now show up as a singularity 
in the residual relative momentum, hence as an infrared effect. 

We also expect that the Lagrangian ([T3D ensures the existence of bound 
states, corresponding to the "unperturbed" heavy quarkonia states. However, 
this also shows that - unlike for heavy-light systems - there is no infinite mass 
limit for quarkonia; the "unperturbed" states described by (|13D will still be 
mass dependent. 

The Lagrangian (|13]) is the minimum that is needed to shift the phases 
from the velocity dependence into the residual momenta. The spin depen- 
dent terms appearing as well in order 1/m do not contribute to the infrared 
behaviour. This is well known from QED, where all infrared contributions 
are independent of the spin of the radiating particle; only the total charge 
plays a role. Spin symmetry thus remains unbroken to leading order, and 
spin symmetry breaking effects may be calculated as perturbations without 
encountering infrared problems. 

The leading order Lagrangian (O) resembles very much the one of non- 



relativistic QCD (NRQCD) as formulated by Caswell and Lepage |T2|. They 



suggest an expansion in f/c, where v is the typical relative velocity of the 
two heavy quarks bound in the quarkonium. Although it does not seem that 
the two approaches are completely equivalent 0, they have many common 
features. 

The binding of the two heavy quarks will generate a small non-perturbative 
scale A, which now in general depends on the heavy mass. Although the 
bottomonium and the charmonium are not Coulombic systems, the case of 
a a/r potential is instructive. Neglecting any running of a, the size of such 
a Coulombic system is RBohr = I /{am), which is large compared to the 
Compton wavelength Xq = 1/m. However, the small scale l/Rsohr depends 
on the mass such that it does not approach a finite limit as m oo, even 
for running a^. 

The Lagrangian (^) is the starting point of an effective theory treatment 
of heavy quarkonium decays. The processes which may be considered in this 
type of approach are decays of heavy quarkonia, in which the two heavy 
quarks annihilate. The annihilation process is governed by a large energy 
scale set by the heavy quark mass m, while the binding of the quarkonium 
introduces a small scale A. The appearance of these disparate mass scales 
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allows for an effective theory treatment, yielding an expansion in powers of 
A/m of the relevant amplitudes of full QCD. 

Heavy Quark Spin symmetry implies that the "unperturbed" heavy quarko- 
nia states fall into degenerate spin symmetry quartets. For a given orbital 
angular momentum I and radial excitation quantum number n, the four 
states (in the spectroscopic notation "^^^^dj) 

\nHt nH.^i nH, nH,+i] (14) 

form such a spin symmetry quartet. An exception are the S waves (£ = 0), 
for which the three polarization directions of the rr'Si and the u^Sq form the 
spin symmetry quartet. The consequences of this symmetry for transitions 
from an excited quarkonium to the ground states have been investigated 
recently |1T3 |. 



In order to exploit the consequences of the spin symmetry for the transi- 
tion matrix elements we shall use the trace formalism. We denote with \Y() 
the spin symmetry quartet consisting of the spin singlet and the spin triplet 
for a given orbital angular momentum £. The coupling of the heavy-quark 
spins may be represented by the matrices 

jj {y\ = { -^+^5 for the spin singlet 
^ I for the spin triplet ' ^ ' 

Using these representations one may then analyse the spin structure of matrix 
elements for processes involving quarkonia. As a simple example we study a 
matrix element like (^) with \A) being a quarkonium state \Y{J, J^, n, I)). In 
the heavy mass limit we have 

(F(J, J„n,/,5)|/iWr/i(;)|0) = a(n,/)Tr {HyT} (16) 

where a{n, I) is independent of the spin coupling of the heavy quarks. A 
simple consequence of spin symmetry is that the basis of the 16 Dirac matrices 
may be reduced to four, which are the generalization of the Pauli matrices 
[H, Furthermore, from (ITO) we have for the ground state quarkonia 



(y(l, J„0,0,l)|/i(+V(</.)/ilr^|0) = (F(0, 0,0, 0,0)1/1^75/^(7)10). (17) 

In a simple wave function picture this means that the wave functions at the 
origin of the two ground state quarkonia are equal in the heavy mass limit. 
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3 Higher Order Correction in HQQET 



In a similar way as in HQET the corrections in HQQET fall into two classes. 
The first type of corrections are the QCD radiative corrections, which in 
general may be calculated systematically. The second class are the recoil 
corrections appearing as a power series in A/m; they originate from matrix 
elements of higher dimensional operators induced by the expansion of the 
currents and the Lagrangian. 

Let us first consider the QCD radiative corrections, which may be calcu- 
lated systematically using the Feynman rules of the effective theory. Based 
on our choice for the leading order term Cq, which is the sum of -^static 
and the first order term h{iD)'^h/{2m), one derives for the propagator of the 
heavy particle with velocity v and residual momentum k 

the corresponding expression for the antiparticle is obtained by the replace- 
ment V — i> —V. This expression contains all orders in 1/m; however, the 
reason why we had to include these higher-order terms was that this removes 
the divergent phase occurring in the limit of small relative velocity. This 
phase may be absorbed as a long-distance effect into the states, which thus 
have to evolve according to the dynamics dictated by Co given in (|13]). The 
"true" short distance contribution is well behaved as the relative velocity 
of the heavy particles vanishes. Furthermore, it may be expanded again in 
powers of (1/m)" without encountering a problem. This is at least true at 
the one-loop level, where (|T8|) is the propagator of NRQCD [|T^, and it has 
been shown in |^ that the coefficients of the ultraviolet divergences at the 
one-loop level may be expanded in 1/m. 

On the other hand, if we start directly from the static limit using the 
propagator of HQET 

H,Uk) = P+^:^, (19) 
vk + le 

imaginary parts will show up, which become ill-defined as velocities of heavy 
quarks coincide, such as 



/ 



'^^^ ^;6{vk)6{v'k) ^ / S^T^iSivk)^ asv-^v' (20) 



(27r)4p V ^ V y J (27r)4fc 



8 



at the one-loop level, which are contributions to the divergent phase. How- 
ever, these are absorbed into the states as being a long distance contribution 
which is generated by the infrared dynamics of Cq and consequently they 
have to be dropped here. 

With this additional prescription, namely to shift the terms diverging 
as the two velocities become equal into the states, we may as well obtain 
a 1/m expansion of the short distance contribution by calculating directly 
with ([T9|), i.e. within the HQET framework. This should coincide with what 
is obtained in using ( |18|) and subsequent expansion in powers of 1/m. The 
diverging terms reappear in the states, since by an appropriate choice of the 
leading order Lagrangian (as in (|13])) the corresponding terms are reproduced 
by the residual momenta, but now they are buried in the non-perturbative 
infrared physics of the states. The price to pay is a mass dependence of the 
states, which is not accessible via a 1/m expansion. 

This simplified prescription allows us to apply the full machinery of HQET 
to the calculation of the short- distance corrections in HQQET, and at the 
one-loop level the results obtained for inclusive heavy quarkonia decays co- 
incide with what is obtained in NRQCD and subsequent expansion in 
powers of 1/m 0. A proof whether this is true also in higher orders in the 
loop expansion lies beyond the scope of the present article. 

The other type of corrections are the power corrections, i.e. the higher 
order corrections in the 1/m expansion. These are induced by the expansion 
of the field and of the Lagrangian and are included as perturbations. It has 
been noticed that the terms of order 1 / and higher of the Lagrangian and 
the fields depend on the convention: As in any effective theory it is possible 
to perform local field redefinitions, thereby shifting certain terms from the 
Lagrangian into the fields and vice versa |T^. However, any physical matrix 



element remains unchanged, if such a field redefinition is performed. For our 
purpose a convenient definition of the higher order terms in the Lagrangian 
is the one, in which all terms, which would vanish by the use of the static 
equation of motion derived from (^, are shifted into the definition of the 
fields. In this convention the Lagrangian up to and including 1/m^ terms 
takes the form 

^ = ^static + (21) 

where 

-^static = ^^^^ ^^^^ - ^^"^ ^^^^ 
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= (^) (A'l + Gi) + f-i-)' (K, + G,) + 0(l/m') , (22) 

\2m/ \2m/ 

Since now only a single velocity f appears, we omit the subscript v from the 
field operators in the following. Furthermore, we have defined 

/,(±)[(,/})2 _ (i^/})2]/,(±) 

m[{tD,),[i-tvD),itD^)]]h^^^ 
{-i)U^^a,,{{iD^), \{-ivD), (.23) 

The corresponding expansion of the field Q^^-* reads 

V 2m 8m^ 

8m^ ^ ^ / 

In fact this is the form that has been obtained from QCD by a sequence of 
Foldy-Wouthuysen transformations in 

The corrections to currents involving heavy quarks are then obtained by 
inserting the expansion (^) for the full fields of QCD, while the corrections 
to the states are implemented via time-ordered products of the expansion of 
the currents with higher-order terms of the Lagrangian. However, here we 
have to take into account the fact that the leading order dynamics of the 
states already contains the first order kinetic energy contribution Ki and 
hence this piece must not to be included as a perturbation. 

The higher order terms of the heavy mass expansion for heavy quarko- 
nia are parametrized then by matrix elements involving the "unperturbed" 
quarkonia states that are described by Cq and operators expressed in terms 
of the static fields hl'^\ Due to the mass dependence of the states fiavour 
symmetry is broken, which means that these matrix elements are different 
for charmonium and bottomonium states. These matrix elements are gen- 
uinely non-perturbative quantities and spin symmetry may be used to count 
the number of independent parameters. Accessing the actual values of these 



Ki-- 






Gi 


= g[^^ + g[-^ 


Gp) = 


K2-- 




4^^ = 


G2 


= + Gir^ 


Gt^ = 
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parameters requires input beyond HQQET; they may be extracted from ex- 
periment, calculated using some model framework, or eventually obtained 
from a lattice calculation. 

4 Conclusions 

The heavy mass limit has been used with great success for systems involving 
a single heavy quark. In this limit the mass gap between the particles and 
the antiparticles becomes infinitely large and, as a consequence, particle and 
antiparticle numbers are separately conserved. HQET has been formulated in 
the one (anti)particle sector, and this is sufficient for almost all applications 
considered so far. 

In order to describe heavy quarkonia one has to deal with the particle- 
antiparticle sector. Treating both heavy constituents in the static limit is only 
possible if the two velocities v and v' are very different, i.e if vv' — 1 = 0{1). 
The naively calculated QCD radiative corrections exhibit logarithmically di- 
vergent and purely imaginary contributions, which behave as 1/ \J {yv'Y — 1- 
Interpreting these pieces in the standard way as contributions to the anoma- 
lous dimension yields then phases, which are ill behaved as f — > v' . 

It has been shown that these phases may be removed by a redefinition of 
the states 0]. The phases are thus a property of the states and and hence 
an infrared effect. After redefinition of the states, the "true" short distance 
contributions, i.e the anomalous dimensions become real. 

In a heavy quarkonium state the velocities of the two heavy constituents 
differ only by a small amount of the order 1/m, and it is appropriate to 
choose the same velocity for both heavy particles, which is then identified 
with the velocity of the heavy quarkonium. However, such a limit may not 
be taken from the expression obtained in the static limit with v ^ v' . 

In the present paper we have shown that one may formulate a Heavy 
Quarkonium Effective Theory (HQQET), despite of these difficulties. The 
key-point is that the leading order term has to be chosen in such a way 
that the diverging phases are generated by its dynamics. This forced us 
to include also the first order kinetic energy term into the leading order 
Lagrangian (0), which determines the evolution of the states. Consequently 
the "unperturbed" states have to depend on the mass and heavy fiavour 
symmetry is lost. In other words, from this point of view it is likely that 
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there is no mass independent, static limit of a heavy quarkonium state, since 
the purely static Lagrangian does not generate the diverging phases and 
probably also does not generate bound states. 

However, the Lagrangian ([T3|) still has an unbroken heavy quark spin 
symmetry, and hence the "unperturbed" heavy quarkonia have to fall into 
spin symmetry quartets, since all four orientations of the two heavy quark 
spins will lead to degenerate states. This symmetry also allows to reduce the 
number of independent non-perturbative parameters. 

This effective theory framework allows us also the calculation of correc- 
tions. The main result of the present analysis is that the short distance 
corrections may be calculated using the methods of HQET. However, if the 
velocities of the two heavy quarks become equal, divergent terms will appear, 
which are infared contributions and which have to be shifted into the states. 
The appropriately chosen Cq (O) contains these divergent terms as piece of 



its infrared dynamics. The remaining expressions are the "true" short dis- 
tance contributions, which remain well behaved as f — > v' . This prescription 
has been applied at the one- loop level and yields the same 1/m expansion as 
a NRQCD calculation with a subsequent 1/m expansion. 

The second type of corrections corresponds to the 1/m expansion of the 
Lagrangian and the fields. Including these contributions proceeds along the 
same lines as in HQET, with the only difference that no time-ordered prod- 
ucts with the first order kinetic energy operator Ki are present, since this is 
included already in the definition of the "unperturbed" states. 

HQQET may be applied to all processes in which the two heavy parti- 
cles annihilate, such as inclusive and exclusive decays of charmonium and 
bottomonium, (electromagnetic and strong) and also to weak decays of B^. 
Some applications of this idea have been investigated in and using related 
methods in ^ , but HQQET opens a wide field of further applications. 
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